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Abstract 

In hot non-Abelian gauge theories, processes characterized by the momentum scale 
g'^T (such as electroweak baryon number violation in the very early universe) are 
non-perturbative. An effective theory for the soft (|p| ~ g'^T) field modes is obtained 
by integrating out momenta larger than g^T. Starting from the hard thermal loop 
effective theory, which is the result of integrating out the scale T, it is shown how to 
integrate out the scale gT in an expansion in the gauge coupling g. At leading order 
in g, one obtains Vlasov-Boltzmann equations for the soft field modes, which contain 
a Gaussian noise and a collision term. The 2-point function of the noise and the 
collision term are explicitly calculated in a leading logarithmic approximation. In 
this approximation the Boltzmann equation is solved. The resulting effective theory 
for the soft field modes is described by a Langevin equation. It determines the 
parametric form of the hot baryon number violation rate as T = ng^^ \og{l / g)T^ , 
and it allows for a calculation of n on the lattice. 
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1 Introduction 



Due to the chiral anomaly, baryon number is not conserved in the stan- 
dard electroweak theory [1]. At zero temperature, the rate for baryon num- 
ber changing processes is exponentially small. It becomes unsuppressed when 
T is of order of the electroweak phase transition or cross-over temperature 
Tc ~ 100 GeV [2]. Knowledge of this rate is a prerequisite for any attempt to 
understand the baryon asymmetry of the universe [3]. 
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Baryon number changing processes are due to topology changing transitions 
of the weak SU(2) gauge fields. In the high temperature phase (T > Tc) these 
transitions are unsuppressed only when the participating gauge fields are 
transverse and have a wavelength of order {g'^T)~^ [4], where g is the weak 
gauge coupling. This is the scale at which finite temperature perturbation 
theory for non-Abelian gauge theories breaks down [5,6]. Therefore, the rate 
in the high temperature phase, the so-called hot sphaleron rate, cannot be 
computed using weak coupling methods. 

Close to thermal equilibrium, the rate for topology changing transitions can 
be written in terms of a different-time correlation function of the form 



C{t,-t,) = {0{h)0{h)), (1.1) 

where (■ ■ ■) denotes the average over a thermal ensemble [7]. The operator 0{t) 
is a gauge invariant function of the gauge fields A^(t, x) at time t. Computing 
such correlation functions non-perturbatively is a difficult task. Euclidean lat- 
tice simulation are of no use because this is a problem with real (Minkowski) 
time. 



In this paper it is shown that correlation functions like (1.1) are determined by 
an effective classical field theory which is described by Eqs. (9.45) and (9.49). 
The great advantage classical field theories have over quantum field theories at 
finite temperature is that one can treat them non-perturbatively, like, e.g. on 
a lattice. The use of classical field theory for computing the hot sphaleron rate 
was suggested more than ten years ago [8]. The "soft" field modes, i.e. modes 
with spatial momenta of order g'^T, are classical because they contain a large 
number of quanta as can be estimated from the Bose distribution function[3] 



However, not all field modes which are relevant to the problem are classical. 
The "hard" {p ~ T) modes have occupation number of order 1 and they 
strongly affect the non-perturbative gauge field dynamics by Landau damping 
which was realized by Arnold, Son, and Yaffe [9]. 

In order to be able to use the classical field approximation, one has to integrate 
out the hard modes which can be done in perturbation theory. At leading or- 
der one obtains the well known hard thermal loops [10]. Time independent 
problems can be described in the imaginary time formalism. Then, the only 
effect of hard thermal loops is the Debye screening of electric interactions on 

^ For spatial vectors I use the notation p = |p| . Four- vectors are denoted by capitals, 
-P^ = {Po, p) and I use the metric P'^ = Pq— p^. 
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the length scale {gT)~^. For time dependent problems hard thermal loops af- 
fect both electric (longitudinal) and magnetic (transverse) degrees of freedom. 
The Landau damping (or dynamical screening) is described by a discontinuity 
of the hard thermal loop resummed propagator for space-like momenta. In 
Ref. [11], it was realized that hard thermal loop induced interactions between 
soft and "semi-hard" {p ~ gT) modes have an effect on the soft dynamics 
which is even stronger than Landau damping by a logarithm of 1/g. These 
interactions determine the characteristic time scale associated with the soft 
non-perturbative dynamics as 



g'T logil/gY ^^-^^ 



The hard thermal loop effective theory describes field modes with p <^ T 
and should therefore be classical. Nevertheless, it is difficult to use for non- 
perturbative lattice calculations. Like any classical field theory at finite tem- 
perature it is plagued with Rayleigh- Jeans UV divergences which are the same 
as in classical finite temperature Yang-Mills theory. For equal-time correlation 
functions the hard thermal loop effective theory reduces to a 3-dimensional 
Euclidean field theory which corresponds to dimensional reduction at lowest 
order. In this case, the divergences can be removed by local counter-terms. 
For different-time correlation functions this is not possible. There are linear 
divergences which are non-local in space and time and on the lattice they even 
break rotational invariance [12,13] .[^ This is in sharp contrast with Euclidean 
lattice theories, where rotational non-invariant terms vanish in the limit of 
small lattice spacing.[^ 

Fortunately, one can also integrate out the momentum scale gT in perturbation 
theory. In this way one obtains an effective theory for the soft field modes only. 
Since the dynamics of the soft modes is strongly damped, this effective theory 
has a much better UV behaviour than classical Yang-Mills theory. In a leading 
logarithmic approximation [11] it is even UV finite [14]. It has recently been 
used for a lattice calculation of the hot sphaleron rate by Moore [15]. 

The most direct approach to integrating out the scale gT is to compute loop 



^ This problem does not arise in scalar (/)^ theory [17,18] where the ultraviolet 
divergences in different time correlation functions are the same as in dimensional 
reduction. 

^ There are proposals for using the hard thermal loop effective theory for lattice 
simulations [12,19]. In [20], the method suggested in [19] was used to compute the 
hot sphaleron rate [20] and the sphaleron rate in the symmetry broken phase [21]. 
One has to use a large numerical value of the Debye mass such that the physical 
hard thermal loops dominate over the UV divergences. This means that one cannot 
take the continuum limit at fixed Debye mass. 
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diagrams within the hard thermal loop effective theory with soft external 
momenta and semi-hard loop momenta. This has been done at one loop level 
in Ref. [16]. However, the one loop approximation is not sufficient for obtaining 
the correct effective theory for the soft field modes. Computing higher loop 
diagrams in the hard thermal loop effective theory is quite tedious. It turns out 
to be much more efficient to use the formulation of hard thermal loop effective 
theory in terms of kinetic equations, which is the subject of this paper. 

The main results of this paper and a brief description of the method employed 
have been presented in Ref. [11] (for a brief review, see [22]). Here, the cal- 
culation is described in detail with emphasis on the accuracy of the various 
approximations. An argument which was used in [11] when solving the (lead- 
ing log) Boltzmann equation (9.27) was not rigorous. This is improved upon 
in Sec. 9.4, while the result of [11] is not affected. An alternative (and also 
rigorous) method for solving (9.27) has been presented by Arnold, Son, and 
Yaffe [23] . In [14] these authors have described a simple physical picture for the 
results of Ref. [11], using the concept of color conductivity [24,25]. Recently, 
the Boltzmann equation was also obtained in [14], [26, 27]. A recent detailed 
derivation of the Boltzmann equation can be found in [28] . 

The paper is organized as follows. The starting point is the hard thermal 
loop effective theory for gauge field modes with momenta small compared 
with the temperature. It can be conveniently formulated in terms of kinetic 
equations derived by Blaizot and lancu and independently by Nair, which are 
briefly reviewed in Sec. 2. In order to integrate out the scale gT, the flelds 
in the kinetic equations are decomposed into soft and semi-hard components 
(Sec. 3). The effect of the semi-hard fields is described by the term ^ in the 
equations of motion for the soft fields, Eq. (3.10). The calculation of ^ is the 
subject of Sects. 4-8. First, the equations of motion for the semi-hard fields 
are formally solved in Sec. 4. There are two types of approximations which 
will be made, these are quahtatively described in Sec. 5. The first type of 
approximation is discussed in Sec. 6 for a scalar toy model. For the second 
approximation one needs to know the characteristic amplitudes of the soft 
fields (Sec. 7). Then, the perturbative expansion of ^ is discussed in Sec. 8. At 
leading order in the equation of motion for the soft fields is described by 
Vlasov-Boltzmann equations containing a Gaussian noise and a collision term. 
In Sec. 9 the 2-point function of the noise and the collision term are explicitly 
calculated in a leading logarithmic approximation. In this approximation the 
Boltzmann equation is solved, and one obtains a Langevin equation for the 
soft gauge fields. It determines the (leading log) parametric form of the hot 
sphaleron rate (Sec. 10). 
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2 Hard thermal loops 



The first step towards an effective classical theory for the soft gauge fields 
is to integrate out the hard field modes. Since the hard modes are weakly 
interacting, this can be done in the one-loop approximation. The dominant 
contributions are obtained when one momentum in the loop is on shell, = 
0.0 For the remaining propagators one can use the large energy (or eikonal) 
approximation 



(g + P)2 2g-P + p2 2Q-P 2qvP' ^ ' ^ 

where P is some linear combination of the external momenta and v = Q/q. 
In this way one obtains the so called hard thermal loops [10]. 

In Abelian theories, there is only a hard thermal loop 2-point function, while 
there are hard thermal n-point functions for any n in non-Abelian theories. 
This point is essential for the calculation in this paper. In [11,16] it was shown 
that at one-loop order in the hard thermal loop effective theory the dominant 
contribution to the effective theory for the soft modes is due to hard thermal 
loop vertices. 

The hard thermal loop effective theory is an effective theory for momenta p 
small compared with T. Therefore, it can be described by classical equations 
of motion. For the present problem this leads to an enormous simplification 
compared with the diagrammatic approach in Ref. [16]. However, the hard 
thermal loop equations of motion [29] are non-local, which is caused by the 
term 1/v-P in the eikonal propagator (2.1). The non- locality is due to the fact 
that the hard particles, which are the quanta of the hard modes, can move 
undisturbed over distances as large as the typical size of the low momentum 
fields. 

Fortunately, there is a local formulation of the hard thermal loop equations of 
motion which was obtained by Blaizot and lancu and by Nair [30]- [33] (see also 
[34]). It is the non-Abelian generalization of the linearized Vlasov equations 
for an electro-magnetic plasma (see e.g. [35]). In addition to the gauge fields, 
these equations contain fields W^lx, v) describing the fiuctuations of the phase 
space density of hard particles; v is the 3- velocity of the particles with = 1. 
The W"' transform under the adjoint representation of the gauge group. The 
first equation is the Maxwell equation for the gauge fields 



We consider only particles with masses much smaller than the temperature. 
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I 



v''W{x,y) 



(2.2) 



where = — igA^ is the covariant derivative. Furthermore, An{x) = 
A1j^{x)T°' with Hermitian generators of the fundamental representation which 
satisfy [T", T*] = if"-bcrpc ^^^^ which are normahzed such that Tr(T°'T'') — 
{1/2)6"-^. All fields couphng to the gauge fields contribute to the hard thermal 
loops in the same functional form. In a SU(A'") gauge theory with Uf chiral 
fermions (or 2n/ Dirac fermions) and Ng scalars in the fundamental represen- 
tation, one has ml = (1/3) (iV + iV,/2 + nf/A)g'^T^. The rhs of (2.2) is the 
current due to the hard particles. The integral in Eq.(2.2) is over the directions 
of the unit vector v. The equation of motion for W reads 



where E is the electric field strength tensor. The conserved Hamiltonian cor- 
responding to Eqs. (2.2)-(2.3) is [32,33] 



Different-time correlation functions like (1.1) can then be computed as follows: 

(1) Compute the solution A{x) of the non-Abelian Vlasov equations (2.2)- 
(2.3) for given initial conditions Ain(x), Ein(x) and l^in(x, v) a,t t — 0. 

(2) Insert the solution into (9[A(tj)]. 

(3) Then C(ti— ^2) is given by the product C[74(ti)]C[y4(t2)] averaged over the 
initial conditions with the Boltzmann weight exp(— /3i7[74in, Ejn, VTm])- 



3 Separating soft and semi-hard degrees of freedom 

The fields in the kinetic equations (2.2) and (2.3) contain Fourier components 
with momenta of order gT and g'^T. The goal of this paper is to obtain equa- 
tions of motion for the soft fields only. We introduce a separation scale /i such 
that 



The fields A, E and W are decomposed into soft and semi-hard modes. 



[v ■ D,W{x,^^)] = v-E(a;), 



(2.3) 




^^T < < gT. 



(3.1) 



6 



A^A + a, E^E + e, W^W + w. (3.2) 

The soft fieldsQ A, E, and W contain the spatial Fourier components with 
p < fi while the semi-hard ones a, e, and w consist of those with k > ^. 
Both W and w describe fluctuations of the hard particle distribution. W (w) 
is the long (short) wavelength part this distribution varying on length scale 
greater (less) than l/fi. A word on notation: below, soft momenta will always 
be denoted by P, while semi-hard momenta will be referred to by K. 

In order to obtain a gauge invariant effective theory for the soft fields, the 
separation has to be done in a way which does not break gauge invariance. In 
this paper the precise form of the effective theory will be calculated explicitly 
only in a leading logarithmic approximation. For this purpose it is sufficient to 
keep only those terms which depend logarithmically on the separation scale. 
For these terms it is irrelevant how the cutoff is realized exactly. 

In the following, we will obtain two sets of kinetic equations, one for the 
soft fields and one for the semi-hard ones. Due to the non-linear terms these 
are coupled. Then, the semi-hard fields will be eliminated and one obtains 
equations of motion containing only the soft fields. 

After the decomposition (3.2), the field strength tensor becomes 

p^u _^ pt.u ^ jMi. _ -^[^M^ a^]^ (3.3) 

where 

//^^ = [D'^,a'^] - [D^a'^]. (3.4) 

The term F^^^ on the rhs of Eq. (3.3) and the covariant derivatives in Eq. (3.4) 
contain only the soft gauge fields. The Ihs of the Maxwell equation (2.2) be- 
comes 

-ig[D,, a'']] - ig[a^, r] - g^[a^, a"]]. (3.5) 

Inspecting one-loop diagrams in the hard thermal loop effective theory [16] one 
can see that in order to integrate out the semi-hard modes only hard thermal 
loop vertices are relevant at leading order and that diagrams containing tree 
level vertices are suppressed by a factor g (see Fig. 1). In the present calculation 

^ For notational simplicity, no new symbols are introduced for the soft modes. From 
now on ^4, E and W will always refer to the soft fields only. 
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Fig. 1. The interactions which are relevant to the leading order effective theory for 
the soft gauge fields. The hard modes dynamically screen the soft and semi-hard 
fields. The relevant hard-semi-hard interaction is small angle scattering which ex- 
changes color charge between hard particles (see Sec. 9). Soft-semi-hard interactions 
in the full quantum field theory can be neglected. Therefore, the relevant interac- 
tions within the hard thermal loop effective theory are due to hard thermal loop 
vertices. 

tree level vertices correspond to interaction terms in Eq. (3.5). Therefore, at 
leading order in g, one can drop all interaction terms in (3.5) which contain 
semi-hard fields, i.e. all but the first two terms. In addition, the covariant 
derivatives in the second term can be replaced by ordinary derivatives. 

The Ihs of Eq. (2.3) becomes 

[v ■D,W]^[v- D,W] + [v ■ D,w] -ig[v ■ a,W] - ig[v ■ a,w]. (3.6) 

The last term on the rhs contains a product of two semi-hard fields. Therefore, 
it contains both semi-hard and soft Fourier components. The latter will turn 
out to be essential for the equations of motion for the soft fields. 

Assuming that the scales gT and g'^T are well separated, the decomposition 
(3.2) is invariant under the infinitesimal gauge transformation 



when the gauge parameter uj^x) has only soft Fourier components. Under (3.7), 
the other fields transform in the adjoint representation. 



(3.7) 
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SW — —ig[W,u!], Sa — —ig[a,uj], 5w — —ig[w,cu]. 



(3.8) 



3. 1 The equations of motion for the soft fields 



Dropping the self-interactions of the semi-hard fields, the soft momentum part 
of the Maxwell equation (2.2) becomes 



[D„F^^{x)]^ml J ^v''W{x,^r). (3.9) 
The semi-hard fields enter the equation of motion for W, 

[v ■ D, W{x, v)] = V • E{x) + ^{x, v), (3.10) 
where ^ is defined as 

C{x, v) = ig[v ■ a{x),w{x, v)]soh- (3.11) 

The subscript "soft" indicates that only spatial Fourier components with p < n 
are included. Under the gauge transformation (3.7), ^ transforms as the fields 
in Eq. (3.8). Therefore, Eqs. (3.9) and (3.10) are covariant under (3.7). 

3.2 The equations of motion for the semi-hard fields 



Since the semi-hard fields are perturbative, one can neglect their self-inter- 
actions. Then, the equations of motion for the semi-hard fields become linear 
in a and w, 



d'^a" - ff'dX^rnl f ^v^'w, (3.12) 

J 4:7T 



[v ■ D,w]-ig[v ■ a,W]^v -e. (3.13) 

In the next Section these equations will formally be solved. For that purpose 
it is necessary to separate the free parts in (3.13) from interaction terms. We 
rewrite (3.13) as 



vdw^v-e + h, (3.14) 
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where the interaction terms have been collected in 



h = ig^[v ■a,W] + [v ■ A, 



w 



(3.15) 



Eqs. (3.12) and (3.14) will be used to write ^ in terms of the soft fields only. In 
(3.12), (3.14) covariance under the gauge transformation (3.7) is not manifest. 
Only if ^ is consistently expanded in powers of g, gauge covariance will be 
recovered. 



4 "Solving" the equations of motion for the semi-hard fields 

In this Section we will obtain an expansion of the semi-hard fields in powers of 
the soft fields starting from the equations of motion in Sec. 3.2. This expansion 
is only a formal solution to the equations of motion since the soft fields them- 
selves depend on the semi-hard ones. The use of this expansion will become 
apparent in the Sec. 6. Our method will first be illustrated with a scalar toy 
model. The extension to the gauge theory case will be straightforward. We will 
restrict ourselves to the transverse semi-hard fields because the longitudinal 
fields do not contribute in the leading logarithmic approximation in Sec. 9. 

Consider a scalar field which is split into soft and semi-hard Fourier compo- 
nents $ and 0. For the equation of motion for (j), corresponding to (3.12) and 
(3.14), we write 



where V is some linear differential operator acting on x. We will also assume 
that enters the equation of motion for $ in analogy with (3.10). The form 
of the equation of motion for $ is not relevant to the present discussion. 

Eq. (4.1) is equivalent to 



where 0o is a solution to the free equation of motion P0o — 0. The propagator 
A satisfies T>A{x — x') = 5'^{x — x'). Both 0o and the propagator A are fixed 
by imposing initial conditions for 0. We will work in momentum space, where 
Eq. (4.2) reads 



V(t){x) = g(l){x)^{x), 



(4.1) 




(4.2) 




(4.3) 
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where k ~ gT and p ~ g'^T. By iterating Eq. (4.3) one obtains a series 



= 00 + 01 + 02 + ■ ■ ■ , (4.4) 

where 



UK) = gA{K) I ^0n-i(^ - PMP)- (4.5) 

Each 0„ is hnear in 0o. Furthermore, 0^ contains n powers of g^. If the soft 
field $ was a fixed background field, Eq. (4.4) would be a solution to the 
equation of motion. However, this is not the case here. $ itself depends on 
because the latter enters the equation of motion for $. 

The gauge theory case can be treated in full analogy with the above example. 
We consider only the transverse gauge fields. In the leading logarithmic ap- 
proximation of Sec. 9 only these will be required. In spatial momentum space, 
the transverse field is defined as 



alit,k)=Vl\k)a^it,k), (4.6) 
where the transverse projector is 



n\k)=5^^-^. (4.7) 



For notational simplicity, the semi-hard gauge fields will be written without the 
subscript "t" in the following. The equations of motion in spatial momentum 
space are then 



a{t, k) + k'^a{t, k)=ml [ ^Vtw{t, k, v), (4.8) 
where vl = Vl''{k)vK The equation for w reads 

w{t, k, v) + ■ kw{t, k, v) = -V ■ ai{t, k) + h{t, k, v). (4.9) 



The equations of motion can be solved by Laplace transformation.^ The 
Laplace transform of a function f{t) is 



^ Or one sided Fourier transformation, cf. Ref. [35]. 
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/(A;°) = J dte'^''f{t). (4.10) 



It defines an analytic function in the upper half of the complex A;°-plane. Tak- 
ing the Laplace transform of Eqs. (4.8) and (4.9) one obtains (cf. Appendix A) 

ein(k)+iA;oain(k)-XV^)=m2 j ^vt«;(X,v), (4.11) 



- 'w;in(k, \)-iv Kw(K, v) = v • ain(k) + ik\ ■ a(K) + h(K, v), (4.12) 

where the subscript "in" denotes the initial values at t = 0. Eqs. (4.11) and 
(4.12) can be "solved" analogous to Eq. (4.3). For this purpose, it is convenient 
to write them in the form 



-K'a\K)-ml j ^^>(i^, v) = £l(i^ ) , (4.13) 
- ikow ■ Si{K) -iv ■ Kw{K,v)^£2{K,y), (4. 14) 

with 



?Uir) = -eL(k)-^fcoaL(k), (4.15) 



i2{K, v) = wUK v) + V • ain(k) + h{K, v). (4.16) 

First, one solves Eq. (4.14) for w. Inserting the result in (4.13), the latter can 
be solved for a which gives 



a^{K) = A\{{K)ei{K) + J ^AUK,v)i2{K,v). (4.17) 
Then one can determine w, 

w{K,w) = Al,{K,Y)e\{K) + J ^A22{K,v,vMK,Y,). (4.18) 
The propagators in (4.17) and (4.18) are defined as 
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A^A{K)^VHk)A,{K), (4.19) 
A\,{K,^^)=ml-^A,{K)vl, (4.20) 

^ii{K,v)^-—-A,{K)vi (4.21) 

A,.(i^, V, vO =4vr5(-^)(v - v')^ - ml^^^^A,{K). (4.22) 

Here At{K) denotes the hard thermal loop resummed propagator for the trans- 
verse gauge fields 



with the transverse hard thermal loop self-energy 



= (4.24) 
Furthermore, S^^^^ is the delta-function on the two dimensional unit sphere, 

JdQ.,f{V)6^'''\^r-V) = f{v). (4.25) 
Inserting (4.15), (4.16) into Eqs. (4.17), (4.18) we obtain 

a'{K) = al{K) + J ^A\,iK, ^)h{K, v), (4.26) 



w{K, v) = w^{K, v) + j ^A22{K, V, Yi)h{K, vi). (4.27) 

This is the gauge theory analogue of Eq. (4.3). ao and Wq are the solutions to 
the equations of motion for h ~ 0, i.e. without the soft "background" . They 
only depend on the initial values at t — 0, 



<iK) = -^^e {A\{iK) - Ai{(0, k)) aiik) - Aii(X)e4(k) 

+ 1 ^A\,{K,v)w,^{k,v), (4.28) 



Wo{K, v) = -^eAi,{K, v)aL(k) - Ai,{K, v)eL(k) 
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+ J ^^22{K, V, v,)w,niK Vi). (4.29) 

In Sec. 9 we will need the 2-point functions of cIq and wq; the results are listed 
in Appendix B. 

By iterating Eqs. (4.26) and (4.27) one obtains an expansion 



a^ao + ai + a2-\ , (4.30) 

w^wo + wi + W2-\ , (4-31) 

where for n > 1 



<{K)^ J -^A{,iK,v)KiK,v), (4.32) 

Wn{K,v) = I ^A22{K,V,Yi)hn{K,Vi). (4.33) 



The term hn contains the Laplace transform of a product of fields. Using 
Eq. (A. 2) it can be written as 



(27r) 

0<imp°<imfc° 

-vati(i^-P)W^'=(P,v)]. (4.34) 

Each term in (4.34) contains only the free solutions ao and wq together with 
the full A and W . Since the equations of motion (3.12) and (3.13) are linear 
in a and w, the a„ and Wn are linear in ao and Wq. 

Inserting the expansion of a and w into Eq. (3.11), one obtains an expansion 
for e, 



e = Co + 6 + ^2 + ---. (4.35) 

Each term in (4.35) is bilinear in the free fields ao and wq. The term ^„ is 
of n-th order in the soft fields. Furthermore, it is of (n+l)-th order in the 
coupling g. This does does not yet tell anything about the actual magnitude 
of which obviously depends on the amplitudes of the soft fields. 
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5 The approximations for ^ 

In Sec. 3.1 we have seen that the effect of the semi-hard modes can be described 
by the term ^ in the equations of motion for the soft fields. In Sec. 4 we have 
obtained an expansion of ^. In the following, it will be shown that one can 
make two types of approximations for ^ which both correspond to an expansion 
in g. 

(1) The first type of approximation, which will be discussed in Sec. 6, is 
the factorization of the thermal averages over initial conditions for the 
semi-hard fields. The reason why this factorization holds is that the cor- 
relation lengths and times of the semi-hard fields are much smaller than 
the corresponding ones for the soft fields. 

(2) The second approximation which will be made is the truncation of the 
series expansion (4.35). The accuracy of this approximation is controlled 
by the characteristic amplitudes of the soft fields which will be discussed 
in Sec. 7. In Sec. 8 it will be shown that at leading order in g only the 
first two terms in (4.35) need to be kept. 

It should be noted that in order to go beyond leading order in g^ one would 
have to take into account terms which were already dropped in Sec. 3 and 
which are not included in the expansion (4.35). In addition, one would have 
to consider corrections to the hard thermal loop approximation for the hard 
modes. 



6 Factorization 

In order to understand the approximation (i), it is instructive to go back to 
the scalar field example which was discussed in the beginning of Sec. 4. This 
will simplify the notation significantly and the extension to the gauge theory 
case will again be straightforward. In analogy with (3.11) we consider 




(6.1) 



For the thermal average over initial conditions for 0o we write 



{<i>o{K)<i>,{K')) =TK{K){2nY5\K + K'), 



(6.2) 
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instead of (B.9)-(B.11).[^ We will further assume that the characteristic fre- 
quency of is of order gT. 



The expansion of (p in Eq. (4.5) yields an expansion of ^ analogous to (4.35). 
The first term in this expansion is 



UP) = 9 J 0j-MK)MP - K). (6.3) 
The next term contains one soft field and one additional power of the coupling, 



Note that in (6.3), (6.4), and all higher ^„ the field 0o always appears in the 
combination 



X = MK-P)MP'-K), (6.5) 

where K is semi-hard while P and P' are soft. After solving the complete equa- 
tions of motion for (j) and $, the solution will contain products of arbitrarily 
many factors of x- order to compute correlation functions like (1.1), one 
has to perform the thermal average over initial conditions. The key simplifica- 
tion which can be made is the following: it turns out that in general, thermal 
averages of products of x's can be approximated by disconnected parts. 



{Xl---Xn) ^ {Xl)---{Xn)- (6.6) 

But since the Xi already appear in the equation of motion for $, it does not 
make a difference whether one first solves the complete equations of motion 
and then uses (6.6), or whether one replaces 



X ^ ix) (6.7) 

already in the equations of motion for $. With the replacement (6.7) one 
eliminates the semi-hard fields from the equations of motion for the soft fields. 
In this way one "integrates out" the semi-hard fields in the present framework. 



^ In the gauge theory case, the delta function of ko + /cg is the "Laplace delta- 
function" in Eq. (B.8). In this Section we will only be concerned about the para- 
metric form of ^, for which this does not play a role. 
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The approximation (6.7) is similar in spirit to the derivation of kinetic equa- 
tions by Khmontovich [37] .[^ The sum of soft and semi-hard fields corresponds 
to "microscopic" and the soft field corresponds to "macroscopic" degrees of 
freedom. The soft fields can be considered as a spatial (not statistical) aver- 
age of the sum of soft and semi-hard fields. Using this analogy, the Boltzmann 
equation (9.27) has been recently obtained in Refs. [26,27]. 

Now let us see why (6.6) holds. Let us first estimate the thermal average of 
one single x. 



{x)=T{2ny6\P-P')A{K-P). 

In the final result for a non-perturbative 
will be integrated over the soft momenta 
eats up one of these integrals. Therefore 



(6.8) 

correlation function of soft fields, (x) 
P and P'. The delta-function in (6.8) 
, we can estimate 



S\P-P')^p,\g'T)-', (6.9) 

where po is the characteristic frequency associated with $. In A{K — P) we 
can neglect the soft momentum p relative to k which is semi-hard. As long 
as Po ^ gT, we have /cq — Po ~ so that the only mass scale appearing in 
A{K — P) is gT. Since A has mass dimension —3, we can estimate A{K — P) ~ 
{gT)~^. In Sec. 9 we will see that the integration over semi-hard momenta gives 
a logarithmic dependence on the separation scale /i. Here we will ignore these 
logarithms, the results of this Section are valid beyond the leading logarithmic 
order considered in Sec. 9. Thus, we have 

(X) ~ Tp^\g'T)-\gT)-' ~ p-o'g-'T-\ (6.10) 

This expression has to be compared with the connected 2-point function of 

two x's. 



{XiX2)conn=T\2Tif5\Pi + P2 - Pi - P!,)Hk^ - P2)A(P{ - K^) 

'5\K^ + K2-P1- P2) + 5\K^ -K2-P1 + P^)] , (6.11) 

which contains only one delta-function of soft momenta only. Proceeding as 
for Eq. (6.10), we estimate 

(XiX2)conn. ~ T%\g'T)-\gT)-'' ~ Po'g-''T-'\ (6.12) 
Therefore, the connected part is suppressed relative to the disconnected part. 



See also §51 of Ref. [35]. 
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(XlX2)conn. ~ {^^J^^ iXl) {X2) ■ (6.13) 

The form of this suppression is easy to understand. The connected 2-point 
function of X1X2 contains only one delta-function for soft momenta, while 
the disconnected part has two. The latter has one additional integration over 
semi-hard momenta / d'^k ~ [qTY- Thus, the suppression factor is 



w ■ 

The same suppression occurs in higher correlation functions 



(Xl ' ' ■ Xn)conn. I fl' j (Xl ' ' ' Xm)conn. (Xm+1 ' ' ■Xn)conn.- (6.14) 

The result of this consideration is that, in general, one can approximate 



in^Un)), (6.15) 

where ((^n)) is defined as ^„ with x replaced by its thermal average. 

There is one qualitative difference between the scalar model discussed in this 
Section and the gauge theory. In the gauge theory case the replacement (6.7) 
gives zero due to the anti-symmetry of the structure constants. Therefore, is 
approximately determined by its connected 2-point function and the discussion 
of .^0 will be postponed to Sec. 8.1. Instead, consider the term ^i. After the 
replacement (6.7) it becomes 



/a k 
j^^A{K)A{K-P). (6.16) 

Assuming that po <^ gT, one can neglect P in the integral which then becomes 
a dimensionless constant {A.{K) and A.{K) have mass dimension —1 and —3, 
respectively). The accuracy of this approximation is of order g. Then ((^i)) 
takes the form 



iUP))) - g'T^P), (6.17) 
or, in coordinate space, 

Uii^))) g'T^x). (6.18) 
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Proceeding similarly for n > 2, one finds that {{Cn{x))) is proportional to 
(<l>(a;))". Furthermore, ((^n)) contains n + 1 powers of g and one factor T due 
to Eq. (6.8). To get the correct dimension, there must be 1 — n powers of gT 
which is the only scale available. Thus {{^n{x))) is of the form 



Obviously, the expansion parameter in the expansion of ^ depends on the 
amphtude of the soft field (or fields in the gauge theory case). This will be the 
subject of the following Section. 



7 The characteristic amplitudes and frequencies of the soft fields 

In this Section we will estimate the typical amplitudes of the soft fields, both 
in coordinate and in momentum space. Since the separation (3.2) was defined 
in terms of spatial momenta, the various soft fields can and do have different 
characteristic frequencies. In Sec. 8 these estimates will be used to determine 
when the series (4.35) can be truncated. Strictly speaking, some estimates can 
be made only a posteriori when the effective theory for the soft dynamics is 
known. Then, one can go back and check the consistency of the approxima- 
tions. Some of these estimates, together with instructive illustrations can be 
found in Ref. [9]. 

Any estimate for the amplitude of A is gauge dependent. However, as long 
as one works with a gauge fixing parameter which does not contain inverse 
powers of the coupling, the amplitude of the soft gauge field at a fixed time 
can be estimated from the tree level propagator of the transverse gauge field 
which is unscreened. This gives 



Non-linear contributions to the propagator cannot be neglected, but they do 
not change the order of magnitude of (7.1). Thus, from (7.1) one concludes 
that 



Consequently, the two terms in the covariant derivative di — igAi are of the 
same order of magnitude, di ~ gAi, which makes perturbation theory for the 
soft gauge fields impossible. In spatial momentum space (7.2) corresponds to 



(6.19) 




(7.1) 



A(i,x)~^r. 



(7.2) 
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Po 


A{P) 


WiP, v) 












g-6T-3 






g-5j^-3 



Table 1 

Order of magnitude estimates for the amplitudes of the Laplace transformed soft 
fields. 



A{t,p)^g-'T-'. 



(7.3) 



The amplitude of the field W can be estimated from Eq. (B.2), where now the 
momenta are of order g^T. This gives 



{W{t,p,v)W{t,p',v'))r^g-'T-\ 
so that 



(7.4) 



W{t,p,v)r^g-'T-^ 
Thus, we find that in position space 



(7.5) 



W{t,^,Y)r^g^T. 



(7.6) 



Now we consider the frequency spectrum of the soft fields. When po ~ gT, 
one can use perturbation theory to compute A{P) and W{P) from their initial 
values. Consequently, A{P) and W{P) can be estimated from the solution to 
the linearized kinetic equations (see Sec. 4). Forpo ~ o'^T, perturbation theory 
breaks down. Still, the solutions to the linearized kinetic equations give the 
correct order of magnitude estimates. The results are listed in Table 1. 

Inspecting Table 1, one can sec that the difference SA{t, p) = A{t, p) — ^4(0, p), 
which can be estimated as poA{P), is small compared with A(0, p), as long as 
Po ^ 9^T- III order to obtain large values of 5A{t, p), one has to consider smaller 
frequencies. In Ref. [9] the frequency scale associated with large SA{t, p) ~ 
A(0, p) was estimated as 



Po^ g T, 



(7.7) 
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again by considering the linearized kinetic equations. Later on, we will see 
that there is a logarithmic correction to this estimate. In this Section, we are 
only concerned about powers of the coupling and we can use (7.7). The gauge 
fields are saturated at these small frequencies. Therefore we have 



A{P)r.g-^T-' (po~/T). (7.8) 

The small frequency part of W cannot be estimated from the linearized kinetic 
equations. Therefore, one has to jump ahead and use the result from the 
effective theory for the soft field modes (see Sec. 9.6). One finds 



W{P,v) ^g-'T-^ (Po~/T). (7.9) 

Note that this is smaller than A{P) by two powers of the gauge coupling. Prom 
Table 1, one can see that the frequency spectra of A and W are qualitatively 
different. The main contribution to A(t) ~ / dpQA{po) comes from po ~ g^T, 
while W{t) receives its main contributions from po of order gT and g^T. 



8 The perturbative expansion of ^ 



Compared with the scalar model discussed in Sec. 6, the gauge theory case is 
complicated by the following circumstances. 

(1) There are several soft fields that can play the role of $ (see Eq. (3.15)). 

(2) The various soft fields have different characteristic frequencies. The gauge 
field is "slow" in the sense that A{t) is saturated by frequencies of or- 
der g^T. In contrast, W{t) is "fast" since W{t) is mainly determined by 
frequencies of order gT and g'^T. If the equations of motion were linear, 
this would not be a problem. Then, only the low frequency tail of W 
would be relevant to the dynamics of the gauge fields. However, due to 
the non-linearities, the high frequency parts of W might sneak into the 
low frequency parts of the equations of motion. 

(3) For the term the replacement (6.7) gives zero due to the anti-symmetry 
of the structure constants of the gauge group. Therefore, one has to keep 
the term in the equations of motion for the soft fields. Due to the 
factorization property (6.6), correlators of can be approximated by 2- 
point functions. In other words, acts on the soft fields as a Gaussian 
noise. 



21 



8.1 The term 



We are now ready to discuss the first term in the expansion (4.35), 

eo(^, v) = -gr'^ {y ■ 4{x)wl{x, v))^^^^ . (8.1) 
In momentum space, it reads (cf. Eq. (A. 2)) 

Co{P,^) = -gr'" J (04^ • «o(^)^S(^ - ^, v). (8.2) 

0<imfc°<imp° 

In Sec. 6 we have argued that one can approximate the product aoWo by its 
thermal expectation value. However, this replacement gives zero for ,^0) because 
(ooWq) is proportional to 5*"^ which is contracted with the anti-symmetric struc- 
ture constant /"^'^. Therefore, we have to leave in Eq. (3.10) as it stands and 
the solution to the equations of motion for the soft fields will contain products 
of ^o's- However, the thermal averages of these products can be approximated 
by a product of 2-point functions, 

vi)C(^2, V2)) ■ ■ • (eo"-^(a;„_i, v„_i)eo""(x„, v„)) 
-|- permutations(l, . . . , n), (8.3) 

which means that acts as a Gaussian noise. 

For the following it will be important to know the characteristic amplitude 
of ^0- Since (^o) vanishes and because the leading order contribution is deter- 
mined by its 2-point function, the amplitude of is determined by 



Co- (6^0) i (8.4) 
Inspecting Eq. (8.1) one sees that ^o{P) is of the form 

UP) - g J d'kx ~ gigTYx, (8.5) 

where x is a product of semi-hard fields analogous to (6.5). Therefore, one can 
apply the arguments of Sec. 6, 



{UP)UP')) - g\gTf{xx)co... ~ g\gTf ( ) (x)^ (8.6) 
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where we have used Eq. (6.13) in the second step. The characteristic amphtude 
of (x) was determined in Eq. (6.10). Inserting the result gives 



{UP)UP'))-Po'9~'T-', (8.7) 
so that we finally obtain 



1 

-3\ 2 



For the characteristic frequency of the soft non-perturbative dynamics (7.7) 
this gives 



UP)-9-''T-' (po~/T). (8.9) 



We will now determine the expansion parameter which controls the accuracy 
of the Gaussian approximation (8.3) for ^q. By using (8.3) one neglects the 
contribution due to the 3-point function 



[■^ola-point function ~ (666) ^ • (8- 10) 

For the rhs can one can use the arguments of Sec. 6, 



ixf. (8.11) 

Comparing this to Eq. (8.8), we find 

[eo(i^)]3-point function ^ (/f ) ^ UP)- (8-12) 

For the frequency (7.7) this gives 

[6(^^)13 -point function -^9UP) {P^^9'T). (8.13) 

Thus, we find that the Gaussian approximation (8.3) for neglects terms 
which are suppressed by one power of g. 
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8.2 The term 



Now consider the second term in the expansion (4.35) of ^, 

^i(a;, v) =ig{[v- ai{x),wo{x,v)] + [v ■ ao{x),wi{x,v)]}^^^. (8.14) 

The parametric form of ((^i)) was determined in Sec. 6, the result was Eq. (6.17). 
Now we have two soft fields, W and A, which can play the role of $. A non- 
zero contribution due to A would spoil the gauge covariance of the equation of 
motion for the soft fields (3.10). However, as we have argued in Sec. 3, ^ trans- 
forms covariantly under soft gauge transformations and such a non-covariant 
term cannot occur (for an illustration of how these contributions drop out, see 
Appendix C). Thus, at leading order in g we have 

((6)) ~ g'TW. (8.15) 

Taking into account the dependence on v, the most general form linear in W 
is 

((ei(x, v))) = g'T I ^I{v,v^)W{x,v,). (8.16) 
The amplitude of ((i^i)) for the frequency scale po ~ g'^T is of order (cf. Tab. 1) 

{{UP,^)))-9-'T-'. (8.17) 

Thus, we find that ^o{P, v) and {{^i{P, v))) are of the same order of magnitude 
when Po is the characteristic frequency of the soft gauge fields (7.7). 

In Sec. 6 it was argued that the leading order contribution to is given by 

6 ^ m, (8.18) 

while higher correlation functions can be neglected. Here the situation is 
slightly more complicated because we have several soft fields and these have 
different characteristic amplitudes. In particular, the leading order contribu- 
tion to ((^i)) found above is relatively small because W is smaller than A. 
In addition, the soft fields have different characteristic frequencies. Thus, the 
frequency in the estimate (6.9) may be different from po ~ g'^T. 

Still, one can use the arguments of Sec. 6 to show that corrections to (8.18) 
are negligible. Since A{P) is larger than W{P) for any po (see Tab. 1), one 
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can simply ignore W and apply the arguments of Sec. 6 to the field A only. 
This simplifies the discussion because A is dominated by small frequencies of 
order g^T. Imagine that A gives a non- vanishing contribution to ((^^i)). This 
would be of order 



((6(P,v)))^~5-^T-^ (8.19) 

The largest possible contribution (gauge covariant or not) to the connected 
2-point function of can then be estimated using (6.13), 



{Ui))a - (g'^) m))Af ~ g~'T-\ (8.20) 

which implies that 

UP)k,.int function -g-'T-'. (8.21) 

Comparing this with Eq. (8.17) we find 

UnUointiunotion-giiUP)))- (8-22) 

Thus, by using the approximation (8.18) one neglects contributions to which 
are suppressed by at least a factor g. 

8.3 Cnforn>2 

In this Section it is shown that does not contribute to the equations of 
motion for the soft fields at leading order in g when n>2. The estimates made 
in this Section are not meant to be exact. Instead, they are to be understood 
as upper limits. 

Let us first consider the term ((^n))- The largest term which is allowed by 
gauge covariance is 

{{L)) - g'^igTf-^W^. (8.23) 
In momentum space it is of the form (ignoring the dependence on v) 

((e„(P))) ~^"(^r)'-"/ d'p, • • •/ dV-i 

W{Pi) ■ ■ ■ W{Pn-i)W{P - Pi • • • - P„_i). (8.24) 
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Note that not only the small frequency {po ~ g'^T) piece of W contributes to 
this integraL We can estimate d'^piW{Pi) ~ W{x) ~ g'^T. An upper hmit for 
the last factor in (8.24) is W{po ~ g'^T, p) ~ g~^T~^. Consequently, 



{(UP))) - g'^'-'T-^ (8.25) 
which shows that ((Cn(-P))) for n > 2 is in fact negligible relative to and 

Estimates for connected correlation functions of ^„ can be obtained in the 
same way as for in Sec. 8.2. First, we estimate the amplitude of the would- 
be contribution of A to ((Cn))- Proceeding as for Vt^, we find 



{{UP))) A - g'^-'T-'. (8.26) 

The connected 2-point function of ^„ can now be estimated using Eq. (6.13). 
For po ~ g'^T one finds 



[en(P)]2.point function ^ (/f ) ^ {{UP))) A ~ g'{{UP)))A. (8.27) 

Consequently, 

[en(P)]2.point function ^^"-'r-', (8-28) 

which is smaller than and .^i by at least a factor g^ when n >2. 



8.4 The Boltzmann equation for the soft fields 



To summarize this Section, we have found that at leading order in the gauge 
coupling, the equation of motion for the soft fields (3.10) can be approximated 
as 



[v ■ D, W{x, v)] ~ V • E{x) + Co(x, v) + Uiix, v))). (8.29) 

Here ,^0 is a Gaussian noise due to Eq. (8.3). Since the expectation value of ,^0 
vanishes, we also had to keep the second term of the expansion (4.35), which 
can be approximated as in (6.15). The terms ^„ with n >2 can be neglected. 
The form of ((^1)) is given by (8.16), so that Eq. (8.29) can be written as 
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[v ■ D, W{x, v)] ^ V ■ E(x) + ^o{x, v) 

+ ^^t/^J(v,viWx,vO. 



(8.30) 



This equation has the form of a hnearized Boltzmann equation for W, which 
contains a Gaussian noise term. The third term on the rhs plays the role of a 
collision term (for a detailed discussion of this equations see Sec. 9.3). 



9 The leading log approximation 



The noise and the collision term in (8.30) contain contributions which are 
logarithmically sensitive to the separation scale /i in Eq. (3.1). These con- 
tributions will be calculated in this Section. They are due to the transverse 
semi-hard gauge fields. Therefore, one can use the expressions for ao and wq 
which were obtained in Sec. 4 and ignore the longitudinal fields. Then, it will 
be argued that, in order to obtain the effective theory for the soft field modes 
at leading logarithmic order, the logarithm can be replaced by log{l/g). In 
this approximation the Boltzmann equation can be solved. This determines 
the current on the rhs of the Maxwell equation (3.9). For the non-perturbative 
dynamics of the soft modes, the kinetic term in (3.9) can be neglected which 
then leads to the Langevin equation (9.49). 



9.1 The 2-point function of 



To compute the 2-point function of C,o, one can use Eq. (8.2) and insert the 
results for the correlators of ao and wq which are listed in Eqs. (B.9)-(B.ll). It 
is simpler, however, to take advantage of translation invariance in time since 
^0 does not depend on the soft background. Then one can proceed as in Eqs. 
(B.6), (B.7) and write 



{CoiPl, Vl)eS(P2, V2)) = -^rri> [{^oiPu Vl)e(P2, V2)) + (1 2) , (9.1) 

Pi ~r P2 

where the subscript "in" denotes the initial value at t = 0. To evaluate the rhs 
of (9.1) one needs the correlations of ag and Wq with the corresponding initial 
values; these are listed in Eqs. (B.3)-(B.5). One obtainsf™] 



In this Section, "~" indicates that these equation are valid with logarithmic 
accuracy. 
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vi)ei^(p2, V2)) « ^ J(Pi, vi, V2)5'^^(27r)353(pi + P2), (9.2) 

where J is the integral 



J(P, V, V) = ^V,/^ I ;^Ai,(X, V) Ai,(P - X, v) (9.3) 

0<imfe°<impO 

+ (Ai^ (i^ - i^) - Ai^ (0, p - k)) A22 (i^, V, v') I . 

After neglecting the soft momentum P, this integral is logarithmically diver- 
gent for A; — > (see below). Since we are integrating out only fields with k > n, 
we will find a contribution proportional to log(l/;u). The integral is convergent 
in the ultraviolet, so that the only other scale in (9.3) is mo ~ gT. Therefore, 
the /i dependent part of (9.3) will be proportional to \og[gT / n). Here, we will 
only compute this logarithmic contribution to (9.3). 

To see how one obtains the logarithmic contribution to (9.3), consider the 
transverse propagator (4.23). For ~ fc, the transverse self-energy (4.24) is 
of order and there is no infrared divergence. Therefore, we must consider 
\k^\ <^k. In this limit, the self-energy (4.24) can be approximated by 



5\1,{K) ~ -ijm^ ^ (|A;o| « k). (9.4) 
4 k 

Then, the hard thermal loop resummed propagator (4.23) becomes 



where 



A k^ 

7k = — r- 9-6 

Thus, for \k^\ <7k the transverse propagator is of order which means 

that it is unscreened. The infrared divergence occurs only in the transverse 
sector. The longitudinal hard thermal loop self-energy stays of order when 
A;o^0. 

With this approximation, assuming |p°|<5'^T, and neglecting p in Afj^(i^ — P) 
the first term in Eq. (9.3) becomes 
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0<imA;'^<imp'' 

m 



2 f d^'k ji 2 1 



0<imfc°<imp° 

The A;*^ integration contour can be closed in the lower half plane picking up the 
poles at A;° = vi • k and k° — — i7k- Note that this is consistent only as long as 
k"^ <S m^. Only then we have 7k <S k, which is required for \k^\ <^ k. Then, 
we can neglect and p in Eq. (9.7) except for the imaginary part of p^. We 
can thus set p° ie. Proceeding similarly for the second term in Eq. (9.3), 
we find 

- <^ (v. ■ <f ^, . 1^'^ ^, . 1, - V . k - t . k - . J ■ '^-^^ 
The first integral in (9.8) is 



/ 



d^k 1 1 2 _ 1 y /• > 1-1/' 
(2^F V • k - i7k ^* " (2^ J J ^yk - i^^ 



Here, the upper Umit for the k integration was chosen as gT which is where the 
low frequency approximation (9.5) ceases to be valid. For the second integral 
we choose vi as the 3-axis. The 1-axis is chosen such that v lies in the 1-3 
plane. As in Eq. (9.9), the integral over the polar angle is saturated at small 
y — cos{9). Then, we have Vt • Vi,t ~ v • Vi. For the integral over the azimuthal 
angle (/? of k one obtains 



J 1 2/ii 

J v'-k — v-k — ie ^ y^l — (v • v')^' 

The remaining integral over the polar angle gives 



I dy- : — - — . 9.11 

yk + ^7k yk - Z7k A;7k 



-1 
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Both in (9.10) and (9.11) the angular integration is saturated at v • k ~ 
v' ■ k ~ 7fc, which means that the approximation (9.5) was consistent. Again, 
the integral over k gives a logarithm log{gT/ij,), so that we finally obtain 



J(P, v,v')«^ log (^^j7(v,v'), (9.12) 



where 



/(V, V) = -5(^^)(v - V) + \ Z^-/^' (9.13) 
The 2-point function of can therefore be written as 



eo(^,v)e^(P',v'))«-2Ar^log (^) /(v,vO 



d""' — '--{27i)H%p + p'). (9.14) 



In configuration space, this corresponds to 



eo"(x,v)eo'(x', v')) «^ -2^^^ log (^1 I{Yy)5^'5\x-x'). (9.15) 



9.2 The collision term ((Ci)) 



In this Section we approximate .^i according to Eq. (6.15), i.e. we replace the 
products of the free semi-hard fields ao and wq by their expectation values. 
We only need to take into account the piece of which contains the field W. 
The contributions due to A vanish after the replacement (6.15) (see Sec. 8.2 
and Appendix C). Then we have 



' ^ ^ J 27r 4 J 4:71 J 27r 4 ^ ^' ^' 



(27r)4 J An J (2? 

0<imfc''<imp'^ 0<imp"<imfc'' 

fVAl.iK, vi)vi • 4{K - P,)wl{P - K, v) 
+ /^^A22(i^, V, vi)v • a^(P - X)vi • 4iK - Pi)] . (9.16) 

Using the results for the 2-point functions of ao and Wo in (B.9) and (B.IO) 
we find 
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V))) ^ iV/r/l^ / l?k/|^_±_(2.,:M3(p _ p,, 

0<imA;°<imp° 0<imp5<imfc° 



+ A22(is:,v,Vi)^;Vi 



mo 



(AlMi^-Pi)-Al^(0,k-pO) 



(9.17) 



The terms in the curly bracket which do not depend on P vanish because for 
them the A;°-integration contour can be moved to +ioo without hitting any 
cut or pole. The remaining terms in the curly bracket do not depend on Pi. 
Therefore, the p^-contour can be closed in the upper half plane which gives a 
contribution due to the pole at — p^, 

0<imp°<imp'' 

Then we obtain 



((ei(P, v))) « Ng'T J ^^J(P, v,vOW^(P, vO, (9.19) 

where J(P, v, vi) is precisely the integral (9.3) we encountered in the previous 
Subsection. With the approximations discussed there, we finally find 

((ei(P, v))) « Ng^og (^^j I ^/(v,vOW^(P, vi). (9.20) 



9.3 The Boltzmann equation 

With the results obtained above we can write Eq. (8.30) as 

[v ■ D, W{x, v)] f=i V • E(a;) + Co{x, v) 

+ AT/Tlog (^^) I ^/(v,vOW^(x,vi). (9.21) 

Eq. (9.21) is a Boltzmann equation for the soft fluctuations of the particle 
distribution W. The term is independent of the fields and therefore acts as 
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a stochastic force. The third term on the rhs of Eq. (9.21) has the form of a 
coUision term. It is due to the interactions with the semi-hard fields. 

Let us see whether Eq. (9.21) is consistent with the Maxwell equation (3.9) 
for the soft fields. This requires that the current 



r{x) =mlj ^v''W{x, v) (9.22) 

is covariantly conserved. Integrating Eq. (9.21) over the direction of v, the 
third term on the rhs drops out due to 



/^/(v,vO = 0, (9.23) 



An 

and one obtains 



D,J'^{x) = mlJ^U^,v). (9.24) 

That is, the current appears not to be conserved. However, we have argued 
that only the 2-point function of should be relevant to the leading order 
behaviour of the soft gauge fields. But the 2-point function of ^0(2^1, Vi) with 
the rhs of Eq. (9.24) vanishes. 



(^eo(xi,vO |^eo(a:,v)^ = 0, (9.25) 
due to Eq. (9.23). Thus, we can replace 

/^6(^,v)c.O, (9.26) 

so that the current is covariantly conserved within the present approximation. 

Due to the term Eq- (9-21) looks as if it was not gauge covariant. However, 
the 2-point function of ,^0) Eq- (9.15), is invariant under gauge transforma- 
tions of ^0- Therefore, a gauge invariant correlation function computed using 
Eq. (9.21) will not depend on the choice of a gauge after the average over ^0 
has been performed. 

The only spatial momentum scales which are left in the problem are /i and g'^T. 
The field modes we are ultimately interested in are the ones which have only 
momenta of order g'^T. The cutoff dependence on the rhs must drop out after 
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solving the equations of motion for the fields with spatial momenta smaller 
than /i. Thus, after the /z-dependence has cancelled, the logarithm in (9.15) 
must turn into \og{gT / (g'^T)) = \og{l/g). Therefore, at leading logarithmic 
order the soft fields satisfy the Boltzmann equation 

[v ■ D, Wix, v)] = V ■ E(x) + ^o{x, v) 

+ iV/T log (1/^7) / ^/(v,vi)W^(a;,vi). (9.27) 

A collision term similar to the one in (9.27), has been obtained previously by 
Selikhov and Gyulassy [24]. It did not contain the second term in (9.13) which, 
as we have seen above, is necessary for the current (9.22) to be conserved. The 
collision term is larger by two powers of the coupling than the collision term in 
an Abelian plasma. This reflects the fact that collisions in Abelian and non- 
Abelian plasmas are qualitatively different, which was realized by Selikhov 
and Gyulassy [24]. The collision term in (9.27) is dominated by small angle 
scattering. In a non- Abelian plasma, even a single small angle scattering event 
can change the color charge of the hard particles, which is what is "seen" by 
the soft fields.[3 



9.4 Solving the Boltzmann equation 

In the previous Subsection we have obtained a Boltzmann equation valid for 
the soft field modes at leading order in log(l/(y') when the momenta in W{x, v) 
are strictly of order g^T. It will now be shown that one can neglect the Ihs of 
Eq. (9.27) up to terms which are suppressed by a factor {\og{l / g))~^ . This is 
an enormous simplification since it turns (9.27) into a relatively simple integral 
equation. It can then be solved for W, which determines the current on the 
rhs of the Maxwell equation (3.9) for the soft gauge fields. 

It is convenient to use a representation of the field W{x, v) in which the 
collision term is diagonal. Since the second term in (9.13) depends only on 
vi • V2, the collision term commutes with rotations of the velocity variable v. 
Thus, the eigenfunctions of the collision term are the spherical harmonics Yi^ 
and its eigenvalues (see Appendix D) only depend on /. 

The I = 1 projection of the Boltzmann equation is obtained by multiplying 
Eq. (9.27) by and integrating over fiv. In Aq = gauge one obtains 



For a recent discussion of earlier approaches to non- Abelian transport equations, 
see Ref. [14]. 
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doW\x) - ^[D\ W{x)] - [D\ W'^{x)] = ^E\x) + ^(a;) (9.28) 

-i-7V/riog(l/5) W\x). 

Here, 

W{x) = j'^^W{x,wl (9.29) 

W\x) = j '^v'W{x, v), (9.30) 

W'^{x) = / ^ (^ V - ^5'^^ W{x, v), (9.31) 
contain the Z = 0, 1, 2 components of W{x, v), respectively. Similarly, 

eo(^) = /^^^eo(:r,v) (9.32) 
is the I — 1 projection of io{x, v). 

It is easy to sec that the first two terms on the Ihs of Eq. (9.28) can be ne- 
glected. The first term contains a time derivative. We are only interested in 
the low frequency part of (9.28), for which this term is very small. It should 
be noted that the first term can be neglected only if one is interested in the 
dynamics of the gauge fields, since Vr(a;, v) is mainly determined by frequen- 
cies larger than g'^T. Now consider the second term in Eq. (9.28). Using the 
GauB law 

[D,,E\x)] = mlW{x), (9.33) 
it can be estimated as 

[D\ W{x)] ~ \{g^TfE\x) ^ g'E\x), (9.34) 
mj-) 

which is small compared with the term E^ on the rhs. 

In the following, it will become clear that the third term on the Ihs of (9.28) 
can be neglected as well. Therefore, (9.28) can be trivially solved for W\ and 
one obtains 
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which is the result of [11]. It should be noted (and will be explained below) 
that the argument which was used in [11] to justify this approximation was 
not rigorous. 

Before discussing this point more carefully, let us first solve Eq. (9.27) for 
the higher / projections of W{x, v) which is simpler. For this purpose it is 
convenient to expand iy(a;, v) and ^0(2;, v) in spherical harmonics, 

W{x, v) = ^ WiUx)YUy^), (9.36) 

Im 

eo(x, v) = ^ eo,/m(x)rz™(v). (9.37) 

Im 

Consider the {l,m) projection of Eq. (9.21) for / > 2, 



+ ^Ng'Tlog{l/g)WiUx). (9.38) 

Here A; are the eigenvalues of the integral kernel (9.13). For / > 2 they are all 
non-zero (see Appendix D). The Ihs is of order g'^TWi±i^rn±i, while the rhs is 
larger by a factor log{l/g). Thus, to leading logarithmic accuracy the solution 
to (9.38) is 



Eq. (9.38) shows that all Wim are of the same order of magnitude when / > 2. 

Now we return to the discussion of Eq. (9.28). In [11] it was argued that the 
term oc \og{l/g)W'^ is large compared with the Ihs because it is logarithmically 
enhanced. This argument is correct for the transverse projections of W\ but 
it is not correct for the longitudinal part. In [23], the longitudinal projection 
Wl was defined as the part of which contributes to [i5*,Vr*]. Wl enters 
the / = projection of (9.27), which is the equation for current conservation 
discussed in Sec. 9.3. In terms of the fields in (9.28) it reads 



doW{x) - [D\ Wl{x)] = 0. (9.40) 

This equation implies that Wl is very small. In particular, the term propor- 
tional to log{l/g)Wl in the longitudinal projection of (9.28) is much smaller 



I thank Peter Arnold, Dam Son and Larry Yaffe for clarifying discussions about 
this point. 
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than the term [D^ , W'^^]i^ on the Ihs, contrarily to the claim in [11]. However, 
due to (9.39), even the term [D^ , W^^] is smaller than by a logarithm. This 
means that one can in fact neglect the Ihs of Eq. (9.28), which is the asser- 
tion made above and which leads to Eq. (9.35). An alternative derivation of 
Eq. (9.35) was recently presented by Arnold, Son, and Yaffe [23]. 

One may wonder what the longitudinal projection of Eq. (9.28) means. Due to 
the smallness of W^^, the leading terms in the longitudinal projection of (9.27) 
are and ^ql- Therefore, these two terms must cancel each other which gives 



[D\K]^-?,[D\Q. (9.41) 

On first sight this result appears a little puzzling. The purpose of this Section 
was to solve the Boltzmann equation which should determine W{x^ v) as a 
function of E and .^o- Now we have obtained a condition for E itself. The 
reason is that we have made use of the Gaufi law to argue that W{x) can 
be neglected in Eq. (9.28). Therefore, Eqs. (9.35), (9.39), and (9.41) must be 
considered as an approximate solution to the set of equations consisting of 
both the Boltzmann equation and of the GauB law. 



9.5 The Langevin equation 

Inserting the current J* = rn^W^ into the Maxwell equation (3.9) one obtains 

A' + [Dj, F^'] ^ -^A' + C, (9.42) 
where the damping coefficient (or color conductivity, see [14]) is given by 



^ 3 Ng^T\og{l/gy ^ " ^ 

The Gaussian white noise Q is proportional to ^q, 



The 2-point function of C can be obtained from Eqs. (9.15) and (D.3), 

(C{xi)C'\x2)) = 2Tj5'^5''''5\xi - X2). (9.45) 
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The equation of motion (9.42) describes an over-damped system. To see this, 
let us estimate the second term on the Ihs. It contains two covariant derivatives 
of the gauge fields. Each derivative is of order g'^T . Thus, this term can be 
estimated asP^ 

[Dj, F^'] ~ D^A ~ g^T'^A. (9.46) 
For the damping term on the rhs we have 



7A ~ r^TTT' (9-4^) 

log(l/^) t 

where t is the characteristic time scale associated with the soft gauge fields. 
Comparing (9.46) and (9.47), we find 



t-'^g^\og{l/g)T. (9.48) 

Therefore, the second time derivative in (9.42) is negligible and the dynamics 
of the soft modes is correctly described by the Langevin equation 



[Dj,F^'] = --fA' + C- (9.49) 

The approximation of dropping the term A has an important effect. The 
effective theory described by Eqs. (9.45) and (9.49) is no longer sensitive to 
an UV cutoff [14]. In other words, this theory is UV finite (aside from the 
"vacuum energy" which does not affect correlation functions like (1.1)). It 
is equivalent to the stochastic quantization [39,40] of Yang-Mills theory in 3 
Euclidean dimensions (in contrast to the usual stochastic quantization, the 
time variable in (9.49) is the physical time). Zinn- Justin and Zwanziger [38] 
have shown that both equal-time and different-time correlation functions are 
renormalized by the same counter-terms as the static theory, except for a 
possible renormalization of the damping coefficient 7. But 3-dimensional Yang- 
Mills theory is UV finite and even the coefficient 7 is not renormalized in 3+1 
dimensions. 

It is remarkable that within the present approximation the effect of the high 
momentum {k ^ g'^T) modes can be simply described by a local damping 
term and a Gaussian white noise. The hard modes affect the soft dynam- 
ics by Landau damping which, in general, is a non-local effect. The modes 
with g'^T log{l/g) ^k<gT are responsible for interactions between the hard 



This estimate does not rely on perturbation theory. For the soft modes both terms 
in the covariant derivative di — gAi are of the same order because of Ai ~ gT (see 
Sec. 7). 
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Fig. 2. The soft gauge fields correspond to extended field configurations (shaded) 
with a typical size R of order {g^T)~^ . Since the soft fields are changing in time, 
they are associated with a long wavelength electric field. The hard modes act like 
almost free particles moving on light-like trajectories (thick lines). They absorb 
energy from the long wavelength electric field which leads to the damping of the soft 
dynamics. The semi-hard field modes (thin lines) are responsible for color changing 
small angle scattering of the hard particles. In the leading log approximation, the 
typical distance between these scattering events is small compared with R, so that 
the damping becomes effectively local. 

particles such that their mean free path is of order {g'^T \og{l / g))"^ . In the 
limit of very weak coupling, this mean free path is much smaller than the size 
R ~ {g'^T)~^ of a soft gauge field configuration. Consequently, the soft modes 
are damped on a relatively small length scale which means that the damp- 
ing is effectively local (see Fig. 2). This effect is characteristic for non-Abelian 
theories, even though the mean free path (or lifetime) of hard particles has the 
same parametric form in Abelian and non-Abelian theories [41]. However, the 
mean free path is determined by small angle scattering which hardly changes 
the momenta of the hard particles. It is the exchange of color charge which 
makes the damping local (see also the discussion at the end of Sec. 9.3). 



9. 6 The characteristic amplitudes and frequencies of the soft fields revisited 

Now that we know the effective theory for the soft modes, at least at leading 
logarithmic order, we can complete the estimates made in Sec. 7. As in Sec. 7, 
we will not be concerned about logarithms. Eq. (9.48) shows that the estimate 
(7.7) is indeed correct. Prom this, one immediately obtains the estimate for the 
amplitude of A{P) in Eq. (7.8). The amplitude of W{P,\) can be estimated 
from the results of Sec. 9.4, 




(9.50) 
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where Eq. (8.9) has been used in the second step. 



10 The hot sphaleron rate 

The rate for electroweak baryon number violation is proportional to the num- 
ber of topology changing transitions per unit time and unit volume [3]. These 
transitions are unsuppressed only if the participating gauge fields have a typ- 
ical spatial momentum of order g'^T [4]. Furthermore, they require that the 
gauge fields change by an amount as large as the field itself, 5 A ~ A. This is 
the dynamics which is described by Eqs. (9.45) and (9.49). 

Since this effective theory has no dependence on the UV cutoff [14], there is 
only one length scale i? ~ ((yf^T)"^, and only one time scale t ~ {g^ \og{l / g)T)~^ 
left in the problem. Consequently, the hot sphaleron rate can be estimated asf^ 



Tr^^^^g'Hog{l/g)T\ (10.1) 
Therefore, at leading logarithmic order, the hot sphaleron rate has the form 



T = Kg^Hog{l/g)T\ (10.2) 

where k is a non-perturbative coefficient which does not depend on the gauge 
coupling and which has been determined by solving (9.49) on the lattice by 
Moore [15]. 
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Again, this estimate does not rely on perturbation theory. 
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Appendix A 



For this paper to be self-contained, this Appendix contains some basic prop- 
erties of the Laplace transformation (4.10) which are used in the main text. 
The original function f{t) is recovered by the inverse transformation 



m=J^e-^'^f{ko), (Al) 

imfco>0 

where the integration contour goes from — oo -|- i x im/co to oo -|- i x im/co- The 
Laplace transform of a product is 



oo 



/ dte^'°'mg{t) = J ^f{po)9{ko-po). (A.2) 



0<impo<imfeo 

To Laplace transform equations of motion one needs the relations 

oo 

I dte''°'f{t) = -/in - ikofiko), (A.3) 



oo 

I dte^^^'m = -/in + ikofin - klfiko), (A.4) 



where the subscript "in" denotes the initial values a,t t — 0. 
Appendix B 



In the calculation of ((^i)) and of the 2-point function of the noise in Sec. 9 
one encounters correlation functions of the transverse free semi-hard fields. 
This Appendix describes how these correlation functions arc obtained from the 
solutions to the free equations of motion (4.28), (4.29), and the Hamiltonian 
(2.4). 

The correlation functions of the initial values ain and Win are obtained by 
thermal average using the Hamiltonian (2.4). Since we are only interested in 
lowest order perturbation theory, the Hamiltonian can be linearized. The 2- 
point function of the transverse gauge fields is the 3-dimensional gauge field 
propagator 
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T 



(B.l) 



Evaluating the 2-point function of Win, one has to account for the Gaui3 law. 
It can be ignored for the leading logarithmic approximation in Sec. 9, where 
only the transverse semi-hard fields contribute. Then, one can simply use 



«(ki,vi)w|;,(k2,V2))=5«^^(27r)353(ki + k2)47r<5(^')(vi - V2). (B.2) 



m 



D 



Using Eqs. (4.28) and (4.29) together with (B.l) and (B.2) one obtains 



(a-(Xi)af,^(k2)) = -z^[Aii(ii') - Ai^(0,k)]5»''(27r)353(k, + k2), (B.3) 



T 



w;o"(i^i,viX(k2,V2)) = ^A22(i^i,vi,V2)5"'(27r)=^53(ki + k2), (B.4) 



D 



(a-(Xi)<Jk2,v)) = K(^i,v)a?^(k2)) = 



T 



— Al2(i^i,v)5"''(27r)V(ki + k2), 



(B.5) 



where the propagators Aq,^ are given by Eqs. (4.19)-(4.22). 



To compute correlation functions of ao and Wq, it is convenient to use trans- 
lation invariance in time to write 



{al^{xM\x2)) = e{h - h)K{ti - t2,xi)ai'(x2)) 

+ e(t2-ti)(a^^(xi)a^o'(^2-ti,X2)). 

Taking the Laplace transform, this becomes 



(B.6) 



a^,\K,)ai\K,)) 



K{K,)at{k,)) + {d^{\.M\K,)) . (B.7) 



1.0 I 1.0 

For Laplace transformed functions the factor 



(B.8) 
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is the analogue of the deha- function 27!'S{ki+k2) in Fourier space (cf. Eq. (9.18)). 
Combining Eqs. (B.3) and (B.7) we find 



fv]^ ~r 1^2 

^ (A^AiKi) - Ai^(o, kO) + ^ {^Am - AAio, k,)) 



(B.9) 



Proceeding similarly for wq one obtains 




(B.IO) 



and 



(«;o"(i^i,vi)«;^(ir2,V2)) 




(B.ll) 



Appendix C 

In Sec. 8 it was argued that terms proportional to A" do not occur in {{Cn)), 
because they would violate gauge covariance. However, if one writes down the 
perturbative expansion (4.35) for i^, one encounters terms containing A. It is 
interesting to see how these terms cancel after the replacement (6.7). 

Consider h in (3.15) and keep the term A. The corresponding contribution to 
^1 is given by (8.14), where ai and wi are now linear in A. Using (4.32), (4.33) 
one obtains 



Now we insert the expressions for slq and Wq which were obtained in Sec. 4 and 
replace the products of the free semi-hard fields by their thermal expectation 
values according to (6.7). Using Eqs. (B.IO), (B.ll) we obtain 




0<imfe°<impO 0<imp5<imfe' 



[f'^%'A\,{K, w,)wl{K - Pi, y{)wI{P - K, v) 

+ /^^A22(i^, V, vi)v • a^(P - K)w'oiK - Pi, Vi)1 . (C.l) 
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ml J {2nY J An J {2tiY p^-p\ 

0<imfc''<imp" 0<impJ<imA;'' 
■ A\P^)[v'^^^{K, Vi) [A22{K - Pi, V, Vi) + A22(P - i^, V, Vi)] 

- i^22{K, V, vi)^;^ [Ai2(P - X, vi) + Ai2(X - Pi, vi)] }. (C.2) 



As in Eq. (9.17), the terms in the curly bracket which do not depend on P 
vanish and one obtains 



{mp.^)))A 



Ng^T f d^k f dQ 



ml 7(2- 

0<imfe°<imp° 

v'A\^{K, vi) A22(P - K, V, vi) - A22(i^, V, Vi)^;^Ai2(P - K, vi) 



(C.3) 



The two terms in the square bracket cancel after the substitution K ^ P — K 
in the second term. Thus, we find that indeed 



(C.4) 



A similar cancellation (I do not write out the straightforward calculation ex- 
plicitly, instead I am restricting myself to a qualitative description) takes place 
to the 74-dependent terms in 



e2"(P,v)=^/ 



abc 



d^k 



(2nY 

0<imA;° <imp" 

+ V • 3i\{K)wl{P -K,y) + v a^{K)wl{P - K, v) 



(C.5) 



Making repeated use of (4.32) and (4.33), one can write a„ and Wn for n = 1, 2 
in terms of ao, Wq, and A. After the replacement (6.7) none of the three terms 
in (C.5) vanishes individually. Each of them consists of two parts. The first 
part of the first term vanishes because the integration contour can be moved to 
infinity. The second part cancels the first part of the second term. The second 
part of the second term cancels the first part of the third term. Finally, in the 
second part of the third term the integration contour can be moved to infinity, 
so that it does not contribute. 
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Appendix D 



In this Appendix we determine the eigenvalues of the coUision term in (9.27). 
The / = 1 eigenvalue enters the damping coefficient in the Langevin equation 
(9.49). Furthermore, consistency of the solution (9.35) to the / = 1 Boltzmann 
equation (9.28) requires the I > 2 eigenvalues to be non-zero (see Eq. (9.39)). 

Since the integral kernel of the collision term is proportional to /(v, v') in 
Eq. (9.13), we consider the eigenvalue equation 



Eq. (9.23) implies that Aq = 0, which is necessary for the current (9.22) to 
be conserved. To determine for / > 1 it is sufficient to consider the eigen- 
value equation for Y/o(cos 9), which is proportional to the Legendre polynomial 
Pi{cos6). Thus the eigenvalues of the kernel / are determined by 



It is also sufficient to consider one special value for v, a convenient choice is 
V = z. Then, we simply have Pi{z ■ \) = P/(l) = 1, so that 



For odd / we have A; = —1, because in this case the second term in (D.3) 
vanishes. For even / > 2 the integral in (D.3) is always smaller than 1, due to 
Pi{x) < 1. Thus all A; are non-zero and negative when / > 1. 



References 

[1] G. t'Hooft, Phys. Rev. Lett. 37 (1976) 8. 

[2] V.A. Kuzmin, V.A. Rubakov, M.E. Shaposhnikov Phys. Lett. 155B (1985) 36. 

[3] For reviews of electroweak baryon number violation and electroweak 
baryogenesis, see A.G. Cohen, D.B. Kaplan, A.E. Nelson, Ann. Rev. Nucl. 
Part. Sci. 43 (1993) 27, |hep-ph/93022To|; V.A. R ubakov, M. E. Shaposhnikov, 
Usp. Fis. Nauk 166 (1996) 493, |hep-ph/9603208| . 

[4] P. Arnold, L. McLerran, Phys. Rev. D 36 (1987) 581. 




(D.l) 




(D.2) 




(D.3) 



44 



[5] A. D. Linde, Phys. Lett. 96B (1980) 289. 

[6] D. J. Gross, R. D. Pisarski and L. G. Yaffe, Rev. Mod. Phys. 53 (1981) 43. 

[7] S. Yu. Khlebnikov, M. E. Shaposhnikov, Nucl. Phys. B 308 (1988) 885. 

[8] D.Yu. Grigoriev and V.A. Rubakov, Nucl. Phys. B 299 (1988) 67. 

[9] P. Arnold, D. Son and L.G. Yaffe, Phys. Rev. D 55 (1997) 6264. 

[10] E. Braaten and R. Pisarski, Nucl. Phys. B 337 (1990) 569; Phys. Rev. D 45 
(1992) 1827; J. Prenkel and J.C. Taylor, Nucl. Phys. B 334 (1990) 199; J.C. 
Taylor and S.M.H. Wong, Nucl. Phys. B 346 (1990) 115. 

[11] D. Bodeker, Phys. Lett. B 426 (1998) 351. 

[12] D. Bodeker, L. McLerran and A. Smilga, Phys. Rev. D 52 (1995) 4675. 
[13] P. Arnold, Phys. Rev. D 55 (1997) 7781. 



[14] P. Arnold, D. Son and L.G. Yaffe, Phys. Rev. D 59 (1999) 105020 , fiep- 



ph/9810216 



[15] G.D. Moore, MCGILL-98-28, |hep-ph/9810"3T3. 



[16] D. Bodeker, NBLHE-99-04, |hep-ph/9903478 



[17] G. Aarts and J. Smit, Phys. Lett. B 393 (1997) 395. 

[18] W. Buchmiiller and A. Jakovac, Phys. Lett. B 407 (1997) 39. 

[19] C. Hu, B. Miiller, Phys. Lett. B 409 (1997) 377. 

[20] G.D. Moore, C. Hu, B. Miiller, Phys. Rev. D 58 (1998) 045001 (1997) 377. 

[21] G. D. Moore, Phys. Lett. B 439 (1998) 357; Phys. Rev. D 59 (1999) 014503. 

[22] D. Bodeker, Proceedings of the Workshop on thermal field theories and their 
applications, Regensburg, ed. U. Heinz, hep-ph/981146E . 

[23] P. Arnold, D. Son and L.G. Yaffe, UW/PT 98-11, |hep-ph/990l'304 . 

[24] A. Sehkhov, M. Gyulassy, Phys. Lett. B 316 (1993) 373. 

[25] H. Heiselberg, Phys. Rev. Lett. 72 (1994) 3013. 

[26] D. Litim, C. Manuel, ECM-UB-PF-99-04, CERN-TH-99-29, |hep-ph/9902"430| . 



[27] M.A. Valle Basagoiti, EHU-FT/9905, |hep-ph/9903462 



[28] J. P. Blaizot, E. lancu, Saclay-T99/026, CERN-TH/99-71, pIep^h/990338£ . 
[29] R. Jackiw, V.P. Nair, Phys. Rev. D 48 (1993) 4991. 

[30] J. P. Blaizot, E. lancu, Phys. Rev. Lett. 70 (1993) 3376; Nucl.Phys. B 417 
(1994) 608. 



45 



[31] V.P. Nair, Phys. Rev. D 48 (1993) 3432 . 

[32] J. Blaizot and E. lancu, Nucl. Phys. B 421 (1994) 565 . 

[33] V.P.Nair, Phys. Rev. D 50 (1994) 4201. 

[34] P.P. Kehy, Q. Liu, C. Lucchesi and C. Manuel, Phys. Rev. D 50 (1994) 4209. 

[35] E.M. Lifshitz, L.P. Pitaevskii, Physical Kinetics (Pergamon Press, Oxford 
1981). 

[36] J. P. Blaizot, E. lancu, Nucl.Phys. B 421 (1994) 565. 

[37] Yu. L. Klimontovich, Kinetic Theory of Nonideal Gases & Nonideal Plasmas 
(Pergamon Press, Oxford 1982). 

[38] J. Zinn-Justin, D. Zwanziger, Nucl. Phys. B 295 (1988) 297. 

[39] G. Parisi, Y.-S. Wu, Sci. Sin. 24 (1981) 483. 

[40] For a review, see J. Zinn-Justin, Quantum Field Theory and Critical 
Phenomena (Oxford University Press, Oxford 1996). 

[41] R. Pisarski, Phys. Rev. Lett. 63 (1989) 1129; Phys. Rev. D 47 (1993) 5589. 



46 



